On the heat flux and entropy produced by thermal fluctuations 
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We report an experimental and theoretical analysis of the energy exchanged between two con- 
ductors kept at different temperature and coupled by the electric thermal noise. Experimentally we 
determine, as functions of the temperature difference, the heat flux, the out-of- equilibrium variance 
and a conservation law for the fluctuating entropy, which we justify theoretically. The system is 
ruled by the same equations of two Brownian particles kept at different temperatures and coupled 
by an elastic force. Our results set strong constrains on the energy exchanged between coupled 
nano- systems kept at different temperature. 

The fluctuations of thermodynamics variables play an important role for understanding the out-of-equilibrium 
dynamics of small systems [1, 2], such as Brownian particles [3-6], molecular motors [7] and small devices [8]. The 
statistical properties of work, heat and entropy, have been analyzed, within the context of fluctuation theorem [9] and 
stochastic thermodynamics [1, 2], in several experiments on systems in contact with a single heat bath and driven 
out-of-equilibrium by external forces or fields [3-8]. Instead the important case in which the system is driven out-of- 
equilibrium by a temperature difference and the energy exchanges are produced only by the thermal noise has been 
analyzed only theoretically on model systems [10-14] but never in an experiment because of the intrinsic difficulties 
of dealing with large temperature differences in small systems. 

We report here an experimental and theoretical analysis of the statistical properties of the energy exchanged between 
two conductors kept at different temperature and coupled by the electric thermal noise. This system is inspired by 
the proof developed by Nyquist [15] in order to give a theoretical explanation of the measurements of Johnson [16] on 
the thermal noise voltage in conductors. In his proof, assuming thermal equilibrium between the two conductors, he 
deduces the Nyquist noise spectral density, which at those time, well before Fluctuation Dissipation Theorem (FDT), 
was the second example, after the Einstein relation for Brownian motion, relating the dissipation of a system to 
the amplitude of the thermal noise. In this letter we analyze the consequences of removing the Nyquist's equilibrium 
conditions and we study the statistical properties of the energy exchanged between the two conductors kept at different 
temperature. This system is probably among the simplest examples where recent ideas of stochastic thermodynamics 
can be tested but in spite of its simplicity the explanation of the observations is far from trivial. Experimentally we 
determine , as function of the temperature difference, the heat flux, the out-of- equilibrium variance and a conservation 
law for the fluctuating entropy, which we justify theoretically. This system is very general because is ruled by the 
same equations of two Brownian particles kept at different temperatures and coupled by an elastic force [11, 12]. Thus 
it gives more insight into the properties of the heat flux, produced by mechanical coupling, in the famous Feymann 
ratchet [17, 18] widely studied theoretically [12] but never in an experiment. Our results set strong constrains on the 
energy exchanged between coupled nano-systems kept at different temperature. 




FIG. 1: a) Diagram of the circuit. The resistances Ri and R2 are kept at temperature T\ and T2 = 296i^ respectively. They 
are coupled via the capacitance C. The capacitances C\ and C2 schematize the capacitance of the cables and of the amplifiers 
inputs. The voltages V\ and V2 are amplified by the two low noise amplifiers A\ and A2 [27]. b) The circuit in a) is equivalent 
to two Brownian particles (mi and 1712) moving inside two different heat baths at T\ and T 2 . The two particles are trapped by 
two elastic potentials of stiffness K\ and K2 and coupled by a spring of stiffness K (see text and eqs.3,4) The analogy with the 
Feymann ratchet can be made by assuming as done in ref.[12] that the particle ml has an asymmetric shape and in average 
moves faster in one direction than in the opposite one. 

The system we consider is sketched in fig. la). It is constituted by two resistances R\ and i?2, which are kept at 
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different temperature T\ and T 2 respectively. These temperatures are controlled by thermal baths and T 2 is kept fixed 
at 296K whereas Ti can be set at a value between 296K and 88K using a circulation of liquid nitrogen vapors. In the 
figure, the two resistances have been drawn with their associated thermal noise generators r]i and 772, whose power 
spectral density are given by the Nyquist formula \ fji\ 2 = AksRiTi (see eqs.3,4 and ref.[20]). The coupling capacitance 
C controls the electrical power exchanged between the resistances and as a consequence the energy exchanged between 
the two baths. The quantities C\ and C 2 are the capacitances of the circuits and the cables. Two extremely low 
noise amplifiers A\ and A 2 [27] measure the voltage V\ and V 2 across the resistances Ri and R2 respectively. When 
Ti = T2 the system is in equilibrium and exhibits no net energy flux between the two reservoirs. This is indeed the 
condition imposed by Nyquist to prove his formula, and we use it to calibrate the system and check the values of the 
capacitances applying Fluctuation-Dissipation-Theorem (FDT) to the circuit (see ref.[20] for the calibration). In the 
following we will take C = lOOpF, d = 680pF, C 2 = 420pF and R ± = R 2 = lOMfi, if not differently stated. 
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FIG. 2: The joint probability log 10 P(Vi, V2) measured at T\ — 296K equilibrium (a) and out of equilibrium T\ — 88K(b). The 
color scale is indicated on the colorbar on the right side. 



The important quantity to consider here is the joint probability P(Vi, V2), which is plotted in fig. 2a) at Ti = T 2 
and at fig. 2b) at TI = 88K. The fact that the axis of the ellipses defining the contours lines of P(Vi, V 2 ) are inclined 
with respect to the x and y axis indicates that there is a certain correlation between V\ and V 2 . This correlation, 
produced by the electric coupling, plays a major role in determining the mean heat flux between the two reservoirs, as 
we well discuss in the following. The interesting new features occur of course when T\ ^ T 2 . The questions that we 
address for such a system are: What are the heat flux and the entropy production rate ? How the variance of V\ an 
V 2 are modified because of the heat flux ? What is the role of correlation between V\ and V 2 ? We will see that these 
questions are quite relevant and have no obvious answers because of the statistical nature of the energy transfer. 




FIG. 3: a) The probability P(Qi jT ) measured at T\ — 296K (blue line) equilibrium and T\ — 88K (magenta line) out of 
equilibrium. Notice that the peak of the P(Qi, T ) is centered in zero at equilibrium and shifted towards positive value out of 
equilibrium. The amount of the shift is very small and it is about ksT2 b) The measured mean value of (Qi, T ) is a linear 
function of (T 2 — T\ ) . The red points correspond to measures performed with the values of the capacitance C\ , C2 , C given in 
the text and r = 0.2s. The other symbols and colors pertain to different values of these capacitance and other r. The values 
of (Qi,t) have been rescaled by the parameter dependent theoretical prefactor A, which allows the comparison of different 
experimental configurations. The continuous blu line with slope 1 is the theoretical prediction (see eq. 7) In the inset the values 
of < Qi > (at C = lOOOpT) directly measured using P(Qi) (blue square) are compared with those (red circle) obtained from 
the equality < Qi >= [g\ — cr^ eq )/Ri, as discussed in the text. 

We consider the electric power dissipated in the resistance R m with m = 1, 2 which reads (j m = V m i m where i m is 
the current flowing in the resistance m. The integral of the power over a time r is the total energy Q m , dissipated 
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by the resistance in this time interval, i.e. Q m ,r — J t T im V m dt. All the V m and i m can be measured indeed 
im = ic — ic m where %c = C ^^"^ is the current flowing in the capacitance C and ic m = ^m^f^ the current 
flowing in C m . Thus rearranging the terms one finds that Q mjT = W mjT - AU mjT where W ljT = fi CV^dt, 

W 2 , r = Sl +T CV 2 ^ L dt and AZ7 m , r = (Cm 2 +C) (V m (t + r) 2 - V m (t) 2 ) is the potential energy change of the circuit m 
in the time r. Notice that W m are the terms responsible for the energy exchanges as they couple the fluctuations of 
the two circuits. The quantities W\, T and W 2 , T can be identified as the work performed by the circuit 2 on 1 and 
vice- versa [19, 21, 24], respectively. Thus, the quantity Q\^ T (Q 2yT ) can be interpreted as the heat flowing from the 
reservoir 2 to the reservoir 1 (from 1 to 2), in the time interval r, as an effect of the temperature difference. As the 
two variables V m are fluctuating voltages all the other quantities fluctuate too. In fig. 3a) we show the probability 
density function P(Qi,r) 7 at various temperatures: we see that Qi jT is a strongly fluctuating quantity, whose PDF 
has long exponential tails. 

Notice that although for Ti < T 2 the mean value of Qi^ T is positive, instantaneous negative fluctuations can occur, 
i.e., sometimes the heat flux is reversed. The mean values of the dissipated heats are expected to be linear functions 
of the temperature difference AT = T 2 — Ti, i.e. (Qi, T ) = At AT, where A is a parameter dependent quantity, that 
can be obtained explicitly from eqs. 3 and 4, as discussed in ref.[20]. This is confirmed by our results, as shown in 
fig. 3b). Furthermore, the mean values of the dissipated heat satisfy the equality (Q 2 ) = — (Qi), corresponding to an 
energy conservation principle: the power extracted from the bath 2 is dissipated into the bath 1 because of the electric 
coupling. This mean flow produces a change of the variances of V m with respect to the equilibrium value , 
i.e. the Nyquist value measured at the same T m but at AT — 0. Specifically we find cr^ — ( ^rn,eq'^~ ^ Qm ^ -^l ( see 
ref.[20]) which is an extension to the two temperatures case of the Harada-Sasa relation [28]. This result is shown in 

the inset of fig. 3b) where the values of ^(j™^ directly estimated from the experimental data (using the steady state 

P(Qm)) are compared with those obtained from the difference of the variances of V\ measured in equilibrium and 
out-of-equilibrium. The values are comparable within error bars and show that the out-of-equilibrium variances are 
modified just by the heat flux. 




FIG. 4: a) The probability P(AS r ) (dashed lines) and P(AStot) (continuous lines) measured at T\ — 296K (blue line) 
equilibrium and T\ — 88K (green lines) out of equilibrium. Notice that both distributions are centered in zero at equilibrium 
and shifted towards positive value in the out-of-equilibrium. b) (exp(— AS to t)) as a function of T\ at two different r = 0.5s 
and r = 0.1s. c) Symmetry function Sym(AStot) = log [P (AStot )/P(— AStot)] as a function of AStot- The black straight line of 
slope 1 corresponds to the theoretical prediction. 



It is now important to analyze the entropy produced by the total system, circuit plus heat reservoirs. We consider 
first the entropy A5V,t due to the heat exchanged with the reservoirs, which reads A5V,t = Qi,t/Ti + Q 2yT /T 2 . 
This entropy is a fluctuating quantity as both Qi and Q 2 fluctuate, and its average in a time r is (AS r , T ) = 
(Q r ,r) 0-/Ti — I/T2) = At(T 2 — Ti) 2 /(T 2 Ti). However the reservoir entropy AS ryT is not the only component of the 
total entropy production: one has to take into account the entropy variation of the system, due to its dynamical 
evolution. Indeed, the state variables V m also fluctuate as an effect of the thermal noise, and thus, if one measures 
their values at regular time interval, one obtains a "trajectory" in the phase space (Vi(t), V 2 (t)). Thus, following 
Seifert [22], who developed this concept for a single heat bath, one can introduce a trajectory entropy for the evolving 
system S s (t) = —Jzb log P(Vi(£), V 2 (t)) : which extends to non-equilibrium systems the standard Gibbs entropy concept. 
Therefore, when evaluating the total entropy production, one has to take into account the contribution 



AS* 



-k B log 



P(y 1 (t + r),F 2 ft + r)) 



(1) 



It is worth noting that the system we consider is in a non-equilibrium steady state, with a constant external driving 
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AT. Therefore the probability distribution P(Vi, V2) (as shown in fig. 2b)) does not depend explicitly on the time, and 
AS s ^ r is non vanishing whenever the final point of the trajectory is different from the initial one: {V\ (t+r), V^(t+r)) ^ 
(Vi(t), V2(t)). Thus the total entropy reads AS tot, r = AS ryT + AS S)T , where we omit the explicit dependence on t, as 
the system is in a steady-state as discussed above. This entropy has several interesting features. The first one is that 
(AS s , r ) = 0, as a consequence (AS tot) — (AS r ) which grows with increasing AT. The second and most interesting 
result is that independently of AT and of r, the following equality always holds: 

(exp(-AS tot /k B )) = 1, (2) 

for which we find both experimental evidence, as discussed in the following, and provide a theoretical proof, as 
discussed in detail in ref. [20]. Indeed, eq. (2) can be proved by relating the reservoir entropy variation ASV,t to the 
system single trajectory probability. Equation (2) represents an extension to two temperature sources of the result 
obtained for a system in a single heat bath driven out-of-equilibrium by a time dependent mechanical force [6, 22] 
and is the the first experimental result obtained on a system driven by a temperature difference. Eq. (2) implies that 
(AS tot) > 0, as prescribed by the second principle. The other consequence of eq. (2) is that AS tot takes into account 
only the out-of-equilibrium fluctuations. Indeed its PDF is a delta function centered in zero when T\ = T2. 

The measured probabilities P(AS r ) and P(AStot) are shown in fig. 4a). We see that P(AS r ) and P (AS tot) 
are quite different and that the latter is close to a Gaussian and reduces to a Dirac S function in equilibrium, i.e. 
Ti = T2 = 296K (notice that the small broadening observed in equilibrium is just due to unavoidable experimental 
noise and discretization of the experimental PDF). Thus we find here a quantity that is rigorously zero in equilibrium, 
both in average and fluctuations, and its mean value and variance provide a measure of the entropy production. The 
experimental measurements satisfy eq.2 as it is shown in fig. 4b). It is worth to note that eq.2 implies that P (AS tot) 
should satisfy a fluctuation theorem of the form log[P(AS to t) / P(—AS to t)] = AStot/ks, Vr, AT, as discussed exten- 
sively in reference [1, 23]. We clearly see in fig. 4c) that this relation holds for different values of the temperature 
gradient. Thus this experiment clearly establishes a relationship between the mean and the variance of the entropy 
production rate in a system driven out-of-equilibrium by the temperature difference between two thermal baths cou- 
pled by electric noise. Because of the formal analogy with Brownian motion the results also apply to mechanical 
coupling as discussed in the following. 

For the experimental analysis we have not used any theoretical input and we will now give a theoretical interpretation 
of the experimental observations. This will allow us to show the analogy of our system with two interacting Brownian 
particles coupled to two different temperatures, see fig. 1-b). Let q m (m = 1,2) be the charges that have flowed 
through the resistances R mi so the instantaneous current flowing through them is i m = q m . A circuit analysis shows 
that the equations for the charges are: 

C C 

R1Q1 = -fty + te-ft)|+^i (3) 

c c 

#242 = -qi y + (21 ~ <te)x + m ( 4 ) 

where X = C2 C\ + C (C\ + C2) and 7? m is the usual white noise: (viVj) — 2SijkBTiRj5(t — t'). The relationships 
between the measured voltages and the charges are : 

qi = (V 1 -V 2 )C + V 1 C 1 (5) 
q2 = (V 1 -V2)C-V 2 C2 (6) 

Eqs. 3 and 4 are the same of those for the two coupled Brownian particles sketched in fig. lb) by considering q m 
the displacement of the particle m, i m its velocity, K m = 1/C m the stiffness of the spring m and K = 1/C the 
coupling spring. With this analogy we see that our definition of Q m corresponds exactly to the work performed by the 
viscous forces and by the bath on the particle m, and it is consistent with the stochastic thermodynamics definition 
[1, 19, 24-26]. 

Thus our theoretical analysis and the experimental results apply to both interacting mechanical and electrical 
systems coupled to baths at different temperatures. Starting from eqs. (3)- (4), we can prove (see ref. [20]) that eq.2 
is an exact result and that the average dissipated heat rate is 

C 2 AT 

= A(T 2 -T 1 )= L ^, (7) 
with Y = [(C± + C)R\ + (C 2 + C)R 2 ] and A = C 2 /(XY) is the parameter used to rescale the data in fig. 3b). 
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To conclude we have studied experimentally the statistical properties of the energy exchanged between two heat 
baths at different temperature which are coupled by electric thermal noise. We have measured the heat flux, the 
entropy production rate and we have shown the existence of a conservation law for entropy which imposes the 
existence of a fluctuation theorem which is not asymptotic. Our results, which are theoretically proved, are very 
general since the electric system considered here is ruled by the same equations as for two Brownian particles, kept 
at different temperatures and mechanically coupled. Therefore these results set precise constrains on the energy 
exchanged between coupled nano and micro-systems kept at different temperature. We finally mention that for the 
quantity Wi an asymptotic FT can be proved both experimentally and theoretically, and this will be the subject of a 
paper in preparation. 
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I. EXPERIMENTAL DETAILS 



A. Experimental set up 



The electric systems and amplifiers are inside a Faraday cage and mounted on a floating optical table to reduce mechanical 
and acoustical noise. The resistance Ri, which is cooled by liquid Nitrogen vapors, changes of less than 0.1% in the whole 
temperature range. Its temperature is measured by a PT1000 which is inside the same shield of R\. The si gnal V\ and V 2 are 
amplified by two custom designed JFET amplifiers [1] with an input current of lpA and a noise of 0.7nVVHz at frequencies 
larger than 1Hz and increases at 8nV\/Hz at 0.1Hz. The resistances Rl and R2 have been used as input resistances of the 
amplifiers. The two signals V\ and V 2 are amplified 10 4 times and the amplifiers outputs are filtered (at AKHz to avoid aliasing) 
and acquired at 8KHz by a 24 bits by an ADC. We used different sets of C\ , C 2 and C. The values of CI and C2 are essentially 
set by the input capacitance of the amplifiers and by the cable length 680pF < CI < 780pF and 400pF < C2 < 500pF. 
Instead C has been changed from lOOpF to lOOOpF. The system has always been calibrated in equilibrium at T\ — T 2 = 296K 
using the FDT and estimating the spectrum using the values of the capacitances, see next sections. 



B. Noise spectrum of the amplifiers 

The noise spectrum of the amplifiers A\ and A 2 (Fig.l of the main text), measured with a short circuit at the inputs, is plotted 
in fig. S. la) and compared with the spectrum Spi of V\ at T\ — 88. We see that the useful signal is several order of magnitude 
larger than the amplifiers noise. 




FIG. S.l: a) The power spectra Spi of Vi measured at Ti = 88K (blue line) ( C = 100pF, Ci = 680pF, C 2 
spectrum of the amplifier noise, b) The equilibrium spectra Sp±( red line) and Sp 2 (green line) measured at Ti 
with prediction of eqs.S.l and S.2 in order to check the values of the capacitances (Ci , C2). 



= 430pF) is compared to the 
= T 2 = 296 K are compared 
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C. Check of the calibration 

The equilibrium spectra of V\ and V 2 at Ti = T 2 used for calibration of the capacitances are: 

c .. Ak B T 1 R 1 {\+Lo 2 (C 2 R 1 R 2 + R 2 (C 2 + C) 2 )} 

Sl(uJ > = {l-^X Rl R 2 ) 2 +^Y (S ' 1} 

o , , 4fc B T 2 gx[l + u 2 (C 2 Rl R 2 + g|(C7i + C) 2 )] 

" {l-u 2 X Rl R 2 ) 2 +u 2 Y (S - 2) 

This spectra can be easily obtained by applying FDT to the circuit of fig. 1 in the main text. 

The two computed spectra are compared to the measured ones in fig. S.la). This comparison allows us to check the values 
of the capacitances C\ and C 2 which depend on the cable length. We see that the agreement between the prediction and the 
measured power spectra is excellent and the global error on calibration is of the order of 1%. This corresponds exactly to the 
case discussed by Nyquist in which the two resistances at the same temperature are exchanging energy via an electric circuit (C 
in our case). 



D. The power spectra of Vi and V2 out-of-equilibrium 

When T\ / T 2 the power spectra of V\ and V 2 are: 

q , v = 4fe B Tx gi[l + u 2 {C 2 R x R 2 + R 2 (C 2 + C) 2 )] Ak B {T 2 - T x )uj 2 C 2 R 2 R 2 

1[UJ) (1 - uj 2 X RiR 2 ) 2 + u 2 Y (1 - uj 2 X RiR 2 ) 2 + uj 2 Y { ' } 

q , x ^ B T 2 gi[l + u 2 (C 2 R x R 2 + R\(C X + C) 2 )] AksiT.-T^uj 2 C 2 R 2 R X 

2[UJ) (1-uj 2 X R ± R 2 ) 2 + uj 2 Y + (1 - uj 2 X R ± R 2 ) 2 + uj 2 Y { ] 

These equations have been obtained by Fourier transforming eqs. S.7,S.8, solving for Vi(cj) and V 2 (uo) and computing the 
modula. The integral of eqs. S.3 and S.4 gives the variances eq. S.24 directly computed from the distributions. 
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FIG. S.2: a) The power spectra Spi of Vi and Sp 2 of V 2 measured at Ti = 120K and T 2 = 296K ( C = lOOpF, d = 680pF, C 2 = 430pF) 
are compared with the prediction of eq.S.3 and S.4 (dashed lines) b) The corresponding Probability Density Function P(Vi) of Vi (green line) 
and P(V 2 ) of V 2 (blue line) measured at Ti = 120K and T 2 = 296K. Dotted lines are the out-of-equilibrium PDF, whose variance is 
estimated from the measure of the heat flux (see fig. 3 in the main text) and eq.S.24 in the following section. The continuous red line is the 
equilibrium P(V 2 ) at T\ = 296K and the black continuous line corresponds to the equilibrium P(Vi) at T 2 = 120K. 



E. Measure of the equilibrium variance of Vi as a function of Ti 

This measure is necessary to estimate < Qi > starting from the measurement of the variances as explained in fig. 4 of the 
main text. We first measure 0-^ (TI) at T\ — T 2 = 296K. In this way one gets the proportionality constant ol\ = 0-^ eq (Ti = 
T 2 )/T 2 as in equilibrium the variance must be proportional to T\. Thus when T\ < T 2 one gets <r^ eq (Tl) = ol\ Ti. As 
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explained in the main text < Q\ >= (a\{Tl) — cr^ (Tl)) / R\. In fig. S.2b) we compare the measured PDF of V\ and V 2 
with the equilibrium and the out-of-equilibrium distributions as computed by using the theoretical predictions eq.S.24 for the 
variance. 



II. DYNAMICAL EQUATIONS FOR V m AND Q m . 

We want to describe the dynamical evolution of both the electric and thermodynamic variables for the circuit in fig. 1 in the 
main text with a set of coupled Langevin equations. For this purpose we write the Langevin equations governing the dynamical 
evolutions for the voltages across the circuit: 

(Ci + C)Vi = CV 2 + ^-(771 - Vx) (S.5) 

(C 2 + C)V 2 = CV X + ^(772 - V 2 ) (S.6) 
^2 

where we have substituted eqs. (5)-(6) into eqs. (3)-(4) in the main text. We rearrange these equations in a standard form, and 
obtain 

Vi = fi(VuV 2 ) + viim+<ri2m = fi(VuV2)+€i (S.7) 

V 2 = f 2 (V u V 2 ) + a 2im + a 22 rj 2 = f 2 (V u V 2 ) + £ 2 (S.8) 

where the "forces" acting on the circuits read 

f (v v\ v v C 2 R2Vi+C{R2Vi+RiV2) „ „ 

f l{Vl ,V 2) = ai V 1+ a 2 V 2 = - [C2C + Ci{ c 2+c)]RiR2 , (S-9) 

C\R 1 V 2 + C(R 2 V 1 +R 1 V 2 ) 

/a(Vl ' Va) = 7lVl+7i,V2 = " [C 2 C + C 1 (C a + C7)]fliJfa ' (S - 10> 



the coefficients Gij read 



C 2 + C 



XR1 

C 

RlV\2 — Rl&21 — j£ 

Ci + C 

° 22 - ~xr 2 ~' 

and the noises ^ introduced in eqs. (S.7)-(S.8) are now correlated (CiCj) = 29ijS(t — t'), where 
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i2i(C 2 C + Ci(C2 + C)) 2 


+ 


R 2 (C 2 C + d(C 2 + C)) 2 


Ti(C(C 2 + C)) 




r 2 (c(Ci + C7)) 


fli^C + Ci^ + C)) 2 


+ 


fi2(C 2 C + Ci(C2 + C)) 2 


TxC 2 


+ 


r 2 (Ci + c) 2 



and 0i2 = 02i • We now notice that the rate of the dissipated heat in circuit m reads 



(S.ll) 
(S.12) 
(S.13) 



Q m = V m i m = ^(V m - Vm ) = (C m + C7)Vi - CV m , , (S.14) 

where ra' = 2 if m = 1, and m' — 1 if m = 2. The rightmost equality in eq. (S.14) follows immediately from eqs. (S.5)-(S.6). 
So one has a formalism where both the voltages and the dissipated heats are described as stochastic processes, driven by the 
thermal noises rj m . 
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III. PROBABILITY DISTRIBUTION FUNCTION FOR THE VOLTAGES 

We now study the joint probability distribution function (PDF) P(Vi, V 2 ,t), that the system at time t has a voltage drop 
Vi across the resistor Ri and a voltage drop V 2 across the resistor R 2 . As the time evolution of V\ and V 2 is described by 
the Langevin equations (S.7)-(S.8), it can be proved that the time evolution of P(Vi, V 2l t) is governed by the Fokker-Planck 
equation [2] 

r\ r\ r\ 2 Q2 

d t P{V u V 2 ,Qut) = L P{V u V 2 ,t) = -_(/ 1 P)-_(/ 2 P) + ( ? 11 _P + ^ 2 _P 

+29 ™iJinr p (S - 15 > 

dvldv 2 

We are interested in the long time steady state solution of eq. (S.15), which is time independent P(Vi,V 2 ,t — > oo) = 
Pss(Vi, V2). As the deterministic forces in eqs. (S.7)-(S.8) are linear in the variables V\ and V 2 , such a steady state solution 
reads 

27re -(aV?+bV 1 V 2 +cV?) 

P..(V U V 2 ) = (S.16) 

V-b 2 + Aca 



where the coefficients 



X {CiT 2 Y + C\CR2l\ + T 2 (CiRi + CRi + C 2 R 2 )}} 
2 [y 2 T!r 2 + C 2 RiR 2 (T 1 - T 2 ) 2 } 
XC[(C 2 + C]R^T\ + (Ci + C)J?iT 2 ] 

[y 2 TiT 2 + ctRMTi - r 2 ) 2 ] ' 

X {C 2 T{Y + C{GR{T 2 + Ti(Cifli + Ci? 2 + C 2 i? 2 )]} 



2[F 2 T 1 T 2 + C 2 J R lJ R 2 (T 1 - T 5 



2; 



can be obtained by replacing eq. (S.16) into eq. (S.15), and by imposing the steady state condition dtP = 0. We are furthermore 
interested in the unconstrained steady state probabilities Pi, ss (Vi), and P 2tSS (V 2 ), which are obtained as follows 



where the variances read 



Pi,ss(Vi) = ( dV 2 P ss {Vi,V 2 ) = ^2= (S.17) 

P 2 ,ss(V 2 ) = /dVLP„(Vi,V r 2 ) = -^2= (S.18) 
J ^2na^ 

2 _ T^C + C^Y + (T 2 -T 1 )C 2 R 1 



(S.19) 

■_■ ' - v v 1 ^-'x/^ — - ■ • * — ^ /com 
a 2 ^7 (S.20) 



2 _ T 2 (C + C 1 )Y-(T 2 -T 1 )C 2 R 2 



XY 

IV. PROBABILITY DISTRIBUTION FOR THE DISSIPATED HEAT AND AVERAGE RATE 

We start by noticing that the heat injected from the bath 1 is then dissipated in the bath 2 (and vice-versa), and so we expect the 
probability distribution of Qi and Q 2 to be symmetric. Thus in the following, we will only study the probability distribution of 
Qi. We now proceed by introducing the joint probability distribution function of the variables Vi, V 2 , and Qi, $(Vi, V 2 , Qi,t) 
As each of these three variables evolves according to a Langevin equation, the time evolution of their PDF is described by the 
Fokker-Planck equation [3, 4] 

dMVuV 2 , Ql ,t) = -^ (/l *)-^( /a *) + ^^* + ^* + JMia ^* 



9 - {m (Vi*) + + 2, 12 A + I*} 

d 2 



+V r 1 2 r 22 -- 2 $ (S.21) 
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with 

rn = Mil + M12, 

r 12 = A;i(9i2 + M22, 

r 22 = k\d xx + k\d 22 + 2k x k 2 d X2 , 

(S.22) 

and £4 = (Ci + C), k 2 — —C. It is worth noting that the first part of the right hand side of eq. (S.21) is identical to the rhs of 
eq. (S.15). 

We proceed by proving eq. (7) in the main text, expressing the dissipated heat rate as a function of the system parameters. We 
have 

d t (Qi) t = d t JdV 1 dV 2 dQ 1 Q 1 ^(VuV 2 ,Qut) = J dVidVadQi QiWV u V 2 , Q,t) = -r n + i- (Vf) 
C 2 AT 

= -nr, (S-23) 

where we have replaced the time derivative dt$(Vi,V 2 ,Q,t) with the rhs of eq. (S.21) and used the equality (V 2 ) = af, with 
G\ as given by eq. (S.19). This equation corresponds to the one given in the main text. 

We can now obtain the expressions for the variance of V\ and V 2 , as introduced in the main text. Using eq.(S.23) we can 
express eq. (S.19) and eq. (S.20) in the following way: 

° 2 m = ^m,eq+ < Qm > Rm (S.24) 

where = Tm ( c + C rn f ) [ s ^ e equilibrium value of cr^ at < Q m >= 0. 

V. CONSERVATION LAW 

We now turn our attention to eq. (2), and provide a formal proof for it. In order to do this, we derive a relation between 
the reservoir entropy change AS r ^ T and the system dynamics. For simplicity, in the following we divide the time into small 
intervals At: let us assume that the system (the circuit in our case) is in the state V = (Vi, V 2 ) at time t, and let's denote by 
V' = (Vi + AVi, V 2 + AV 2 ) its state at time t + At. Let V F (V -> V'|V, t) be the probability that the system undergoes 
a transition from V to V 7 provided that its state at time t is V, and let Vr{V' — > VIV 7 , t + At) be the probability of the 
time-reverse transition. We have 

P F (V^V'\V,t) = y > dr ?1 dr ?2 ^(AVi-At(/l(Vi,y 2 ) + a 11 7 ?1 +a 12 r ?2 )) 

x5(AV 2 - At(f 2 (V u V 2 ) + crsir/i + cr 22 77 2 ))pi (r/i)p 2 (r7 2 ), (S.25) 
where 5(x) is the Dirac delta function. Given that the noises are Gaussian distributed, their probability distributions read 



PmiVra) exp 



A * (S.26) 



A7TR m k B T 7 



4Rmk B T 

and expressing the Dirac delta in Fourier space S(x) = l/(27r) J dq exp(iqx), eq. (S.25) becomes 

P F (V^V'\V,t) = ^ J d gi d 92 expK 9l Ay 1 + 92 Ay 2 )] y"nd^e A 4^ (/ ~ + ^^ ,B) -^](S.27) 



( At r 

'{ ±k B T{T 2 L 



= exp <^ - C?fliT 2 (Vi - /1) 2 + ClR 2 T 1 (V 2 - f, 



2) 



+2C(V! -h-V 2 + f 2 )(C 1 R 1 T 2 (V 1 - h) - C 2 R 2 T X (V 2 - f 2 )) 



+ C>(R 2Tl + RiT 2 )(Vi - h-V 2 + f 2 f] } ^ V tf 5 (S - 28> 
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where we have taken AVm/ At ~ V m , and exploited the fact that all the integrals in eq. (S.27) are Gaussian integrals. Similarly, 
for the reverse transition we obtain 

Pr(V ^V\V',t + At) = y'dr ?1 d»72 5(Ay 1 + At(/l(y 1 ',V r 2 , ) + (Tiir ? i+ ( T 12 r ?2 )) 

x6(AV 2 + M(f 2 (V{, V£) + a 2im + a 22 r l2 ))p 1 {i ll )p 2 {r l2 ) (S.29) 
= exp|- 4fc ^ T2 [cfR.T^V, + A) 2 + C\R 2 T X {V 2 + f 2 f 

+2C(V! + h-V 2 - f2)(CiRiT 2 (V 1 + h) - C 2 R 2 T X {V 2 + f 2 )) 
+ C\R 2Tl + R^m +fi-V 2 - h f] } ^/ff- (S-30) 

We now consider the ratio between the probability of the forward and backward trajectories, and by substituting the explicit 
definitions of fi{V u V 2 ) and f 2 (Vi,V 2 ), as given by eqs. (S.9)-(S.10), into eqs. (S.28) and (S.30), we finally obtain 

P F (V^V'\V,t) ( (Ci + C)^ - CV 2 {C 2 +C)V 2 -CV 1 \ ( QiQ 2 \ 
l0S P R (V>^V\V>,t + At)=- At { W + W J =A ^W + wJ' (S31) 

where we have exploited eq. (S.14) in order to obtain the rightmost equality. Thus, by taking a trajectory V — )> V over an 
arbitrary time interval [£, t + r], and by integrating the right hand side of eq. (S.31) over such time interval, we finally obtain 

P F (V^V'\V,t) (Q 1 Q 2 \ AQ 
k B lQ g 577^ ^mr, + , Z\ = Ur + 7^ = AS r , r (S.32) 



•Ph(V'^v|V',* + t) V^i r 5 

We now note that the system is in an out-of-equilibrium steady state characterized by a PDF P ss (Vi, V2), and so, along any 
trajectory connecting two points in the phase space V and V the following equality holds 



exp [AStot/ks] = exp [(A5 r , r + AS s , r) /k B ] 
V F (V^V f \V,t)P ss (V) 



(S.33) 



v R (v ^v\v',t + T)p sa (vy 

where we have exploited eq. (S.32), and the definition of AS S:T in eq. (1) in the main text. Thus we finally obtain 

V F (V V'| V, t)P ss (V) exp [-AS tot /k B ] = P R (V V| V, t + r)P ss (V) (S.34) 

and summing up both sides over all the possible trajectories connecting any two points V, V 7 in the phase space, and exploiting 
the normalization condition of the backward probability, namely 

V l V/ ^ + T)Pea(V) = 1, (S.35) 

V',V 

one obtains eq. (2). It is worth noting that the explicit knowledge of P SS (V) is not required, in order to prove eq. (2). 

Finally, we note that, from a general perspective, eqs. (S.7)-(S.8) correspond to the Langevin equations of a stochastic system, 
whose variables V\ and V2 interact through non-conservative forces, and where the white noise is correlated. Therefore our proof 
of eq. (S.35), and thus of eq. (2) in the main text, holds in general for systems with such characteristics. 
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